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Abstract 

We report on measurements of the mass and total decay width of the W boson 
with the L3 detector at LEP. W-pair events produced in e + e~ interactions between 
161 GeV and 183 GeV centre-of-mass energy are selected in a data sample corre- 
sponding to a total luminosity of 76.7 pb _1 . Combining all final states in W-pair 
production, the mass and total decay width of the W boson are determined to be 
M w = 80.61 ± 0.15 GeV and T w = 1-97 ± 0.38 GeV, respectively. 
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1 Introduction 



For the 1997 data taking period, the centre-of-mass energy, y/s, of the e + e~ collider LEP at 
CERN was increased to 183 GeV. This energy is well above the kinematic threshold of W-boson 
pair production, e + e~ — > W + W~. 

Analysis of W-pair production yields important knowledge about the Standard Model of 
electroweak interactions pQ through the measurements of the mass, Mw, and the total decay 
width, T w , of the W boson [2]. These parameters were initially measured at pp colliders jSE]- 

First direct measurements of in e + e~ collisions were derived from total cross section 
measurements [SHH]j mainly at the kinematic threshold of the reaction e + e~ — > W + W~, y/s = 
161 GeV, where the dependence of the W-pair cross section on the W-boson mass is largest. 
At centre-of-mass energies well above the kinematic threshold, the mass and also the total 
width of the W boson are determined by analysing the invariant mass of the W-boson decay 
products HDHISI- 

In this letter we report on an improved determination of the mass and the total width of 
the W boson. The analysis is based on the data sample collected in the year 1997 at an average 
centre-of-mass energy of 183 GeV, corresponding to an integrated luminosity of 55.5 pb" 1 . 
The invariant mass distributions of 588 W-pair events selected at this energy are analysed to 
determine Mw and Tw The results based on the 1997 data are combined with our previously 
published measurements based on the 1996 data collected at centre-of-mass energies of 161 GeV 
and 172 GeV ■ 

2 Analysis of Four-Fermion Production 

During the 1997 run the L3 detector [T3j collected integrated luminosities of 4.04 pb" 1 , 49.58 pb _1 
and 1.85 pb" 1 at centre-of-mass energies of 181.70 GeV, 182.72 GeV and 183.79 GeV, respec- 
tively, where these centre-of-mass energies are known to ±0.05 GeV ^Hj- These data samples 
are collectively referred to as 183 GeV data in the following. 

The W boson decays into a quark-antiquark pair, such as W~ — > ud or cs, or a lepton- 
antilepton pair, W~ — > £~i>t (£ = e, /i, r); in the following denoted as qq, tv or // in general for 
both W + and W~ decays. Four-fermion final states expected in W-pair production are ivtv{pj\ 
qqiv{pi), and qqqq(j), where (7) indicates the possible presence of radiative photons. 

The following Monte Carlo event generators are used to simulate the signal and back- 
ground reactions: KORALW t& and HERWIG [H] (e+e" -> WW -> ffff(j)); EXCAL- 
IBUR [HI (e+e" -> fffffr)); PYTHIA [TJ] (e+e" -> qq(>y), ZZ( 7 )); KORALZ |2D| (e+e" -> 
//+//" (7), T+T "(7)); BHAGENE3 jH], BHWIDE |22j and TEEGG ^ (e+e" -> e+e" (7)), 
DIAG36 and LEP4F |25] (leptonic two-photon collisions); PHOJET [H] (hadronic two- 
photon collisions). The response of the L3 detector is modelled with the GEANT |2Z| detector 
simulation program which includes the effects of energy loss, multiple scattering and showering 
in the detector material. 

The selections of the four-fermion final states are described in detail in References IoTIdI and 28 
for the data collected at yfs = 161 GeV, 172 GeV and 183 GeV. These analyses reconstruct 
the visible fermions in the final state, i.e., electrons, muons, r jets corresponding to the visible 
t decay products, and hadronic jets corresponding to quarks. In order to select a pure sample 
of qqqq events, the cut of 0.67 on the neural-network output described in the qqqq cross-section 
analysis is applied [2H]- Kinematic constraints as discussed below are then imposed to improve 
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the resolution in the measured fermion energies and angles and to determine those not measured. 
The invariant mass of the W boson is obtained from its decay products. 

The mass and the width of the W boson are determined by comparing samples of Monte 
Carlo events to the data. A reweighting procedure is applied to construct Monte Carlo samples 
corresponding to different mass and width values. Using this method, effects of selection and 
resolution are automatically taken into account. 



3 Event Reconstruction imposing Kinematic Constraints 

The final states qqev, qqfiv and qqqq contain at most one primary unmeasured neutrino. For 
each event a kinematic fit is performed in order to determine energy, Ef, polar angle, Of, and 
azimuthal angle, 0/, for all four fermions, /, in the final state. The kinematic fit adjusts 
the measurements of these quantities for the visible fermions according to their experimental 
resolutions to satisfy the constraints imposed, thus improving their resolution. 

Four-momentum conservation and equal mass of the two W bosons are imposed as con- 
straints, allowing the determination of the momentum vector of the unmeasured neutrino. For 
the energy constraint, the exact centre-of-mass energies as given in the previous section are 
used. For hadronic jets, the velocity (3f = \pf\/Ef of the jet is fixed to its measured value as 
systematic effects cancel in the ratio. For qqev and qq\xv events, this yields a two-constraint 
(2C) kinematic fit, whereas for qqqq events it is a five-constraint (5C) kinematic fit. 

Events with badly reconstructed hadronic jets are rejected by requiring that the probability 
of the kinematic fit exceeds 5%. The kinematic fit mainly improves the energy resolution and 
less the angular resolution. The resolutions in average invariant mass, m- 1TW) typically improve 
by a factor of four for qqev and qq^xv events and a factor of six for qqqq events. 

For qqrv events, the decay products of the leptonically decaying W boson contain at least 
two unmeasured neutrinos in the final state. Therefore only the hadronically decaying W boson 
is used in the invariant mass reconstruction. The energies of the two hadronic jets are rescaled 
by a common factor so that the sum of their energies equals half the centre-of-mass energy, thus 
imposing equal mass of the two W bosons. The rescaling improves the resolution in invariant 
mass by nearly a factor of four. Since invariant masses of W bosons in Ivlv events cannot be 
reconstructed as the decay of both W bosons involves neutrinos, Ivlv events are not used in 
the analysis for W mass and width. 



4 Fitting Method for Mass and Width 



The fitting procedure uses the maximum likelihood method to extract values and errors of the 
W-boson mass Mw, and the total width Tw, denoted as \I/ for short in the following. In fits to 
determine M w only, the Standard Model relation T w = 3G F M^/(2 v / 27r)(l + 2a s /(3vr)) gSj is 
imposed. Otherwise, M w and T w are treated as independent quantities. 

The kinematic fit imposing the equal-mass constraint determines the weighted average of the 
two invariant W masses in an event, m inv , which is considered in the fit for mass and width. The 



total likelihood is the product of the normalised differential cross section, L[ 
for all data events. For a given four- fermion final state i, one has: 



minvi ^0, evaluated 



Lj(m inv , ^) 



1 



BG 



dcTi(m inv , dof G (m inv ) 



dm ir 



+ 



dm ir 



(1) 
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where <Tj and af G are the accepted signal and background cross sections and fi(&) a factor 
calculated such that the sum of accepted background and reweighted accepted signal cross 
section coincides with the measured cross section. This way mass and width are determined 
from the shape of the invariant mass distribution only. The total and differential cross sections 
of the accepted background are independent of the parameters \& of interest. They are taken 
from Monte Carlo simulations. 

The total and differential signal cross sections depend on For values \l/ n t varied during 
the fitting procedure, these cross sections are determined by a reweighting procedure applied 
to Monte Carlo events originally generated with parameter values ^ gen - The event weights Ri 
are given by the ratio: 

2 

i2j(Pl,P2,P3,P4,%,W fl t,W gen ) = -j, (2) 



gen; 



where Aii is the matrix element of the four-fermion final state i. The matrix elements are 
calculated for the generated four-vectors, (px,P2,P3,P4, & 7 ), of the four fermions and any ra- 
diative photons in the event. Since the Monte Carlo sample used for reweighting is based on 
the three Feynman graphs in W-pair production (CC03 [2HH3T] ). the matrix element in the 
denominator is calculated using only CC03 graphs. The matrix element in the numerator is 
based on all tree- level graphs contributing to the four-fermion final state i. The calculation of 
matrix elements is done with the EXCALIBUR |TB] event generator. 

The total accepted signal cross section for a given set of parameters ^ nt is then: 



gen 

0-<(*flt) = -^-E^U^nt^gen), ( 3 ) 

« i 

where af en denotes the cross section corresponding to the total Monte Carlo sample containing 
Nf en events. The sum extends over all Monte Carlo events j accepted by the event selection. 

Based on the sample of reweighted events, two methods are used to obtain the accepted dif- 
ferential signal cross section in reconstructed invariant mass m inv . Both methods take detector 
and selection effects as well as ^/-dependent changes of efficiencies and purities properly into 
account. 

In the box method the accepted differential cross section is determined by averaging 
Monte Carlo events inside a mi nv -bin centred around each data event. The size of the bin 
considered is limited by the requirement of including no more than 1000 Monte Carlo events, 
yielding bin sizes of about ±35 MeV at the peak of the invariant mass distribution. In addition, 
the bin size must not be larger than ±250 MeV around m- irw . 

In the spline method, the continuous function describing the accepted differential cross 
section is obtained by using a cubic spline to smooth the binned distribution of reconstructed 
invariant masses. At the kinematic limit of y/s/2 the value of the spline is fixed to zero, while 
at the lower bound of 65 GeV the value of the spline is fixed to the average over a 2 GeV 
interval. The spline contains 25 knots in total. Four knots are placed at each endpoint with 
the remaining knots placed such that an equal number of Monte Carlo events separates each 
knot. 

Both methods yield identical results within 15% of the statistical error. For the numerical 
results quoted in the following, the spline method is used. 

The fit procedure described above determines the parameters without any bias as long as the 
Monte Carlo describes photon radiation and detector effects such as resolution and acceptance 
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functions correctly. By fitting large Monte Carlo samples, typically a hundred times the data, 
the fitting procedure is tested to high accuracy The fits reproduce well the values of the 
parameters of the large Monte Carlo samples being fitted. Also, the fit results do not depend 
on the values of the parameters \l/ gen of the Monte Carlo sample subject to the reweighting 
procedure. 

The reliability of the errors given by the fit is tested by fitting for each final state several 
hundred small Monte Carlo samples, each the size of the data samples. The width of the 
distribution of the fitted central values agrees well with the mean of the distribution of the 
fitted errors. 

5 Mass and Total Width of the W Boson 

Based on the data collected at 172 GeV and at 183 GeV, the mass of the W boson is determined 
for each of the final states qqeu, qqfiv, qqru and qqqq in separate maximum likelihood fits. For 
mass fits in the qqqq channel, the pairing algorithm to assign jets to W bosons used in the event 
selection [HIEE] is changed. The pairing yielding the highest likelihood in the 5C kinematic fit 
is chosen. The fraction of correct pairings is reduced to 60% for the best combination and it is 
25% for the second best combination. However, the signal-to-background ratio in the relevant 
signal region around m inv w 80 GeV is improved. The loss of correct pairings is recovered by 
including the pairing with the second highest likelihood as an additional distribution. Monte- 
Carlo studies show that the two values for Mw obtained from fitting separately the distributions 
of the best and the second best pairing have a correlation of (—1.3 ± 1.0)%, which is negligible. 

The observed invariant mass distributions together with the fit results for the semileptonic 
final states are shown in Figure The distributions of the first and second pairing in qqqq 
events are shown in Figure El while the distribution summed over all final states and both qqqq 
pairings is shown in Figure H3 Combined results are determined by averaging the results of 
individual channels taking statistical and systematic errors into account. The results of fits for 
Mw are summarised in Table ^ The observed statistical errors agree well with the statistical 
errors expected for the size of the data samples used. The results of fits for M w and T w are 
summarised in Tabled 

6 Systematic Effects 

The systematic errors on the fitted W mass and width are summarised in Tables El and EJ They 
arise from various sources and are divided into systematic errors correlated between final states 
and systematic errors uncorrelated between final states. 

6.1 Correlated Errors 

The beam energy of LEP is known with an accuracy of 25 MeV for the 1997 data and 30 MeV 
for the 1996 data, where 25 MeV of these errors are fully correlated j!5j . The relative error on 
Mw is given by the relative error on the LEP beam energy, while the width is less affected. The 
spread in centre-of-mass energy of about 0.2 GeV adds in quadrature to detector resolution and 
total width of the W boson and is thus negligible. 

Systematic uncertainties due to incomplete simulations of initial-state radiation (ISR) are es- 
timated by comparing the Monte Carlo generators KORALW and EXCALIBUR implementing 



5 



different QED radiation schemes. For final-state radiation (FSR), events with FSR simulation 
are compared to events without any FSR and a third of the difference is taken as a systematic 
error. 

The reconstruction of hadronic jets is examined by studying hadronic 55(7) events collected 
at the Z pole and at 183 GeV. A systematic error for the jet measurement is assigned from 
varying the jet energy scale by 0.2 GeV, smearing the jet energies by 5% and smearing the 
jet positions by 0.5°. Effects due to fragmentation and particle decays are determined by 
comparing signal events simulated using string fragmentation as implemented in the PYTHIA 
Monte Carlo program and cluster fragmentation as implemented in the HERWIG Monte Carlo 
program to simulate the hadronisation process. 

The fitting method itself is tested by fitting to various Monte Carlo samples generated 
with known values for Mw and Tw, varying over a range of ±0.5 GeV. The systematic error 
due to the fitting method includes the effects due to different procedures for reweighting and 
smoothing of the invariant mass distributions and choice of technical parameters such as spline 
parameters, box size and occupancy. 

Limited Monte Carlo statistics introduces a tendency of the method to have a slope of the 
linear function relating fitted mass to generated mass less than one. All Monte Carlo samples, 
approximately one million events, are used in the reweighting procedure to minimise this effect 
when fitting data. Fitting several Monte Carlo samples and using the remaining Monte Carlo 
as reference the non-linearity is found to be negligible. 

6.2 Uncorrelated Errors 

The systematic error due to the size of the signal Monte Carlo sample used for reweighting is 
estimated by dividing it into N parts of equal size, N between 2 and 100, and making N fits to 
the same data sample. The spread of the fit results, divided by the square root of N-l, is found 
to be independent of N and yields the systematic error due to Monte Carlo statistics. 

Selection effects are estimated by varying the cut on the probability of the kinematic fit 
and the interval of reconstructed invariant masses being fitted. Effects due to background 
are determined by varying both the total accepted background cross section within its error as 
evaluated for the cross section measurement as well as the shape of the invariant mass spectrum 
arising from the background. 

For qqqq events, strong final state interactions (FSI) between the hadronic systems of the 
two decaying W bosons due to effects of colour-reconnection (HHHS] or Bose-Einstein correla- 
tions may affect the mass reconstruction. In both cases, possible effects are estimated 
by comparing signal simulations including and excluding the modelling of such effects and 
assigning the mass difference found as systematic error. In case of colour reconnection, two 
models, called superconductor model type I and type II as implemented in PYTHIA 5.7 are 
studied [33], adjusted such that they both yield 35% reconnection probability. In case of Bose- 
Einstein correlations, the simulation of this effect as implemented in PYTHIA 5.7 is used [36J. 

For qqeu and qqfiu events, the reconstruction of the lepton energy and angles also affects the 
invariant mass reconstruction. In analogy to hadronic jets, control samples of P. + P.~( r y) events 
selected at the Z pole are used to cross check the reconstruction of leptons. Energy scales and 
resolutions are varied within their errors and the resulting effect on W mass and width is quoted 
as a systematic error. 
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6.3 Z Mass Reconstruction as Consistency Check 



All aspects of the mass measurement, ranging from detector calibration and jet reconstruction 
to fitting method are checked using e + e~ — > qqj events selected at yfs = 183 GeV. For such 
events, the hard initial-state radiative photon reduces the centre-of-mass energy of the e + e~ 
interaction. The presence of the Z resonance causes the distribution of the invariant mass of 
the jet-jet system to exhibit a peak at the Z mass, as it originates from Z decay, with a shape 
similar to the W mass spectrum. 

A kinematic fit is used to improve the mass resolution, enforcing four-momentum conserva- 
tion in order to improve resolutions in energies and angles of measured photons and of the two 
jets and to determine the energy of one photon or two photons escaping along the beam axis. 
For the extraction of the Z mass from the invariant mass spectrum the same method as for the 
W mass measurement is applied. Monte Carlo events are reweighted according to the ratio: 

fl z (^,<,M*J = , (4) 



gen/ 



using the differential cross-section da/dys' where y s' is the reduced centre-of-mass energy 
after initial-state radiation at Monte Carlo generator level. 

The reconstructed mass spectrum together with the fit result is shown in Figure 0] A total 
of 3351 events are selected in a mass window ranging from 70 GeV to 110 GeV. The fitted 
Z-mass value is Mz = 91.172 ± 0.098 GeV, where the error is statistical. Within this error, 
the fitted Z mass agrees well with our measurement of the Z mass derived from cross section 
measurements at centre-of-mass energies close to the Z pole, M z = 91.195 ± 0.009 GeV |3*7j . 
The good agreement represents an important test of the complete mass analysis method. 



7 Results 

The results on M w determined in the qqev, qq^v, and qqru final states are in good agreement 
with each other, as shown in Table H They are averaged taking statistical and systematic 
errors including correlations into account, and compared to the result on Mw determined in 
the qqqq final state, also shown in Tabled The systematic error on the mass derived from qqqq 
events contains a contribution from possible strong FSI effects. Within the statistical accuracy 
of these measurements there is no significant difference between Myv as determined in qqlv and 
qqqq events: 

AM W = M w (qqqq) - M w (qq£u) = 0.35 ± 0.28 (stat.) ± 0.05 (syst.) GeV . (5) 

For the calculation of the systematic error on the mass difference, the systematic errors due to 
strong FSI are not included. 

Averaging the results on Myv obtained from the qqlv and qqqq event samples, including also 
FSI errors, yields: 

M w = 80.58 ± 0.14 (stat.) ± 0.08 (syst.) GeV . (6) 

The summed mass distribution is shown in Figure El and compared to the expectation based on 
this W-mass value. The good agreement between the data and the reweighted mass spectrum 
is quantified by the \ 2 value of 26 for 30 degrees of freedom which corresponds to a probability 
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of 66%. The mass values obtained in fits which determine both Mw and T\y are the same as 
before within 20 MeV while the error on the mass is unchanged. 

Within the statistical error, the width of the W boson determined in qqqq and qqlv events 
agree as shown in Table El For all final states combined the result is: 

T w = 1.97 ± 0.34 (stat.) ±0.17 (syst.) GeV, (7) 

with a correlation coefficient of +10% between Mw and Tw as shown in Figure El Our result 
on Tw is in good agreement with the indirect measurement at pp colliders, 2.07 ±0.06 GeV [I], 
and measurements at LEP [13,38]. It also agrees well with the Standard Model expectation, 
2.08 GeV, calculated for the current world-average W mass [39]. 

The results on Mw presented here agree well with our result derived from the measurements 
of the total W-pair production cross section, M w = 80.781^ (exp.) ± 0.03 (LEP) GeV 0. 
Combining both results yields: 

M w = 80.61 ±0.15 GeV. (8) 

This direct determination of Mw is in good agreement with the direct determination of Mw at 
pp colliders [3J and at LEP at lower centre-of-mass energies |7 H9TITTHT5] and at 183 GeV [SB] ■ It 
also agrees with our indirect determination of Mw at the Z peak, Mw = 80.22 ± 0.22 GeV |37j . 
testing the Standard Model at the level of its electroweak corrections. 
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Table 1: Number of events used in the analysis and results on the mass of the W boson, Mw, 
combining the data collected at 172 GeV and at 183 GeV. The first error is statistical and the 
second systematic. Also shown is the statistical error expected for the size of the data sample 
analysed. 
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Table 2: Results on the mass of the W boson, Mw, and its total decay width, Tw, combining 
the data collected at 172 GeV and at 183 GeV. The first error is statistical and the second 
systematic. Also shown is the correlation coefficient between M w and Tw- 



14 



Svstematic Errors 
on M\ki fMeVl 


Final State 


qqQiz 


nnnu 




nnnn 


LEP Energy 


Zo 


Zo 


Zo 


Zo 


TOD 

IS-tt 


1 K 
10 


10 


10 


10 


r osx 


1 n 
ru 


1 n 

1U 


1 n 

1U 


i n 

1U 




30 

ou 


30 


30 


5 


Fragmentation and b)ecav 


30 


30 


30 


60 


Fitting Method- 


15 


15 


15 


15 


Total Correlated 


55 


55 


55 


69 


MC Statistics 


20 


20 


50 


10 


Colour Reconnection 








70 


Bose-Einstein Effects 








60 


Selection 


20 


20 


20 


20 


Background 


5 


10 


30 


10 


Lepton Measurement 


15 


15 






Total Uncorrelated 


32 


34 


62 


95 


Total Systematic 


63 


64 


82 


118 



Table 3: Systematic errors in the determination of M\y for the different final states. The contri- 
butions listed in the upper part of the table are treated as correlated when combining different 
final states. The contributions listed in the lower part are treated as uncorrelated between 
channels. Total errors are obtained by adding the individual contributions in quadrature. 
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Table 4: Systematic errors in the determination of Tw in qqiv and qqqq production. The 
contributions listed in the upper part of the table are treated as correlated when combining the 
two final states. The contributions listed in the lower part are treated as uncorrelated between 
channels. Total errors are obtained by adding the individual contributions in quadrature. 
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Figure 1: Distributions of reconstructed invariant mass, mi nv , after applying the kinematic fit 
using the equal-mass constraint for events selected in the 183 GeV data: (a) qqev, (b) qqfiu, (c) 
qqrv, (d) qq£v, combining qqev, qqfiv and qqrv. The solid lines show the result of the fits of 
Mw to the indicated final states. The quoted error combines statistical and systematic errors 
in quadrature. 
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Figure 2: Distributions of reconstructed invariant mass, m inv , after applying the kinematic fit 
using the equal-mass constraint for qqqq events selected in the 183 GeV data: (a) first pairing, 
(b) second pairing. The solid lines show the result of the fit of J% to both pairings. The 
quoted error combines statistical and systematic errors in quadrature. 
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Figure 3: Distribution of reconstructed invariant mass, m mv , after applying the kinematic fit 
using the equal-mass constraint for all W-pair events selected in the 183 GeV data used for the 
mass analysis. For qqqq events, both pairings are included. The solid line shows the result of 
the fit of M w . The quoted error combines statistical and systematic errors in quadrature. 
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Figure 4: Distribution of reconstructed invariant mass, m- mv , after applying the kinematic 
fit for qcpi events with hard initial-state radiation selected at 183 GeV. Shown is the region 
corresponding to the radiative return to the Z. The solid line shows the result of the fit of Mi. 
The quoted error is statistical. 
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Figure 5: Contour curves of 68% and 95% probability in the (Mw, Tw) plane from a fit to the 
combined 172 GeV data and 183 GeV data (statistical errors only). The point represents the 
central values of the fit. The Standard Model dependence of IV on Mw is shown as the line. 
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